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Abstract. I discuss galaxy cluster surveys based on the Sunyaev-Zel'dovich 
effect and their relevance for cosmological studies. The unique aspects 
of cluster selection by this method are emphasized and certain issues of 
surveying are addressed. Finally, I briefly present prospects for upcoming 
surveys. 



1. Introduction 

Galaxy cluster catalogs are useful for cosmological studies for a number of rea- 
sons: cluster abundance, distribution and evolution reflect the structure for- 
mation process, which depends sensitively on certain cosmological parameters, 
such as the matter density J7m (Oukbir & Blanchard 1992, 1997; Bahcall 1998; 
Bartlett 1997; Borgani et al. 1999 & these proceedings; Eke et al. 1998; Henry 
1997; Viana &: Liddle 1999); and cluster internal properties provide clues to 
the nature of dark matter and on certain mechanisms associated with galaxy 
formation, such as heating and feedback. 

The construction of well-defined catalogs appropriate for statistical stud- 
ies demands good understanding of the cluster selection criteria (see Postman, 
these proceedings). These criteria in turn depend on the type of observations - 
optical, X-ray, radio - as well as on the data analysis proceedures. Each type of 
observation has its own advantages and disadvantages that we may refer to as 
its intrinsic biases. For example. X-ray cluster selection eliminates the problem 
of projection effects suffered by optical observations, but on the other hand is 
significantly influenced by the rather poorly understood physics taking place in 
the core region. 

The question that arises is, what kind of selection is best? The answer 
clearly depends on the utlimate use of the catalog, and for cosmological studies 
we may argue that the most ideal selection would be based on cluster mass. 
This is because cluster mass and redshift are the two primary theoretical cluster 
descriptors. From this line of argument, understanding cluster selection criteria 
is reduced to relating the observational selection criteria to cluster mass and red- 
shift. This is a useful viewpoint, for it provides a framework for quantifying the 
biases associated with each type of observation for the purposes of cosmological 
studies. We do not argue that one particular observational method is the best; 
rather, that the important point is to use several methods with different biases. 

The Sunyaev-Zel'dovich (SZ) effect (Sunyaev &: Zel'dovich 1972) will soon 
open the way to a new method of surveying for clusters. Several ground-based 
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instruments optimized for SZ surveying are currently under construction, and 
the Planck satellite will supply an almost full-sky catalog of several 10^ clusters 
at the turn of the decade. In anticipation and in the spirit just outlined, we 
shall now examine in more detail this new kind of survey (the 'SZ band'), and 
in particular the nature of its cluster section. 

This potential of surveying for clusters with the SZ effect has been appreci- 
ated for some time. To my knowledge, the first to give source count predictions 
were Korolyov, Sunyaev, & Yakubstev (1986), and many authors have since con- 
sidered the issue. An indicative, useful, though by no means exhaustive, list of 
references focused on SZ surveying (i.e, counts and redshift distributions of SZ 
detected clusters) would include: Bond & Meyers (1991), Bartlett & Silk (1994), 
Markevitch et al. (1994), Barbosa et al. (1996), Eke, Cole, & Frenk (1996), 
Aghanim et al. (1997), Colafrancesco et al. (1997), da Silva et al. (2000), 
Holder et al. (2000), Lo et al. (2000), Benson, Reichardt, & Kamionkowski 
(2001), Fan & Chuieh (2001), Kay, Liddle, & Thomas (2001), Kneissl et al. 
(2001), Springel, White, & Hernquist (2001), Xue & Wu (2001)0. 



2. SZ Cluster Selection 

In the non-relativistic limit (which applies in practice to all but the hottest 
clusters), the SZ effect may be written as 

y = — Trnaxdl 

Jlos 

where Sii, is the brightness change induced in the cosmic microwave background 
(CMB) by compton scattering off electrons (mass m) in the hot intracluster 
medium (ICM) at temperature T and electron density n. The Compton y- 
parameter measures the amplitude of the energy transfer to the photons (ctt is 
the Thompson cross section); and the function ji, describes the spectral form of 
the distortion. It is convenient to consider the effect integrated over the entire 
cluster face: 

M{T) 



Y = dn y{n) oc fg^^s-j^ — (1) 

where now M is the cluster (virial) mass, /gas is the ICM gas mass fraction, (T) 
is the (true) mean electron temperature (i.e., particle weighted), and -Dang(-z) 
is the angular-size distance. Notice that Y is measured in, for instance, square 
arcmins. 

It is from these equations that we begin to understand the nature of SZ 
cluster selection (or detection). Let us emphasize four aspects: 

1. The spectral signature ji, is unique - negative at low frequency, positive 
at high frequency, with a fixed zero point at A ~ 1.4 mm - and the same 



^The reader will appreciate the precarious position of the author in giving such a list. In offering 
my apologies to those that I have unjustly, but without malice, forgotten, I ask no more than 
to be corrected by a kind email in reminder. 
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for all clusters (jj, depends only on the observation frequency). There is 

therefore no k-correction; 

2. The integrated effect oc D~fg, which is the statement that the surface 
brightness is independent of redshift z. This is in contrast to other emission 
mechanisms, that vary as the luminosity distance Aum = (1 + z)^£)ang, 
and which suffer from 'cosmic dimming' (the redshift); 

3. The integrated signal is directly proportional to the total ICM thermal 
energy fgg^M{T) (this is essentially the definition of the mean temperature 
(T)). We may therefore expect a tight correlation between the measurable, 
Y and the cluster mass, e.g., Y oc fgasM^/^, especially in the more massive 
systems where non-graviational effects appear less important. 

4. Cluster detection and subsequent modeling may be based on the same 
measurable quantity Y. 

A particularily Tiscful quantity for characterizing the selection properties of 
a survey is the mininiimi detectable mass as a function of redshift, M^eti^)- For 
example, it determines the source counts as 



where dn/dM is the mass function. A comparison of -M^et from all-sky X-ray 
and SZ (the Planck curve) surveys is made in Figure 1. The X-ray survey flux 
limit corresponds to that of the RASS. An interferometer-like similar to those 
proposed (see below) and capable of surveying several square degrees from the 
ground is also shown. To make the figure, a non-evolving — T relation of 
the form given by Arnaud & Evrard (1999) and a T — M relation from Evrard, 
Metzler, &; Navarro (1996) have been adopted. The cosmological model is an 
open, low-density model, and Eq. (1) was used with a constant gas mass fraction 
/gas = O.OGh-^-^ (Evrard 1997) 

The figure is illustrative, especially because we have no idea of the evolution 
of the Lx^T relation to such rcdshifts. It is nevertheless instructive to note the 
different behavior of the X-ray and SZ curves, and in particular the almost flat 
asymptote of the SZ curves; in fact, they demonstrate a slow decrease towards 
high redshift (this is independent, of course, of the Lx — T relation). In other 
words, SZ selection corresponds essentially to a mass selection. The conclusion 
depends in detail on the constancy of /gas with redshift, but the potential con- 
sequences are interesting: for example, this implies that there is no Fmax for a 
SZ detected cluster. Even more significant is that SZ selection therefore finds 
the same kinds of objects (mass-wise) over a large range of redshifts. For evolu- 
tionary studies of clusters, this is an important advantage over other selection 
techniques that force one to compare ever more massive objects detected at high 
redshift to less massive local samples. 

3. Aspects of Surveying 

The discussion so far as assumed that clusters arc unresolved - that the inte- 
grated SZ fiux of Eq. (1) is directly measured. This situation applies in practice 
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Figure 1. Minimal detectable mass as a function of redshift for a 
Planck-like all-sky survey, an interferometer-like (limited sky, ground- 
based) survey and for a flux-limited X-ray survey similar to the RASS. 
For the SZ curves, the dashed lines account for cluster extension, while 
the solid lines assume that clusters are point sources. The Planck- 
like survey is based on the Planck 2.1 mm channel performance only: 
Y = 10~^ arcmin^/beam and an angular resolution of ^fwhm = 5 ar- 
cmins (this is most surely an optimistic estimate as the Planck SZ 
maps will be weighted by channels with differing resolution and sen- 
sitivity, and foreground removal will further degrade the sensitivity). 
The interferometer-like curve assumes ^fwhm ~ 2 arcmin and a sensi- 
tivity of y = 5 X 10~^ arcmin^. 
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to the Planck survey, with angular resolution of at best 5 armins. Ground-based 
instruments, on the other hand, will operate at higher resolution and clusters 
will appear extended, complicating the selection criteria: some clusters may be 
resolved out, relative to a hypothetical point-source approximation. Detection 
depends now not only on Y, but also on the detailed SZ cluster profile; and 
photometry becomes an issue, as with any extended source. 

Many ground-based instruments will operate, at least at first, at a single 
frequency. In this case, primary CMB fiuctuations (in the background) become 
an important source of noise, and one must extract extended sources over a 
fluctuating background. Another potentially important effect is confusion by 
point sources, either randomly distributed on the sky, or preferentially in cluster 
cores. Without spectral information, their removal requires higher resolution 
imaging; for this reason, some ground-based interferometer projects plan to 
control this confusion with additional telescopes placed on long baselines. 

There are certainly other issues, for example other foregrounds, that we 
could list as problematic. However, I would instead like to address a specific 
question that often arises when preparing a survey: In a given amount of ob- 
serving time, should one go (integrate) 'deep', or instead 'wide', covering more 
sky to less depth? Actually, the real approach is to try to do several surveys, 
i.e., in some sense to do both. All the same, it is useful to try to answer the 
question in terms of maximizing the number of detected objects. This depends 
on the slope of the source counts as a function of survey sensititivit]]^ 

3.1. Deep or Wide? 

We calculate the counts using Eq. (2). For unresolved sources, the detection 
mass is only a function of sensitivity, l^oise, and redshift - Mdet{Ynoise, z). Ad- 
ditional factors are important in the case of resolved source detection, such 
as angular resolution (fipix, the solid angle of a 'pixel') and source profile: 
Mdet{Ynoise, z,Clpix, pro file). A simple analytic model of extended source de- 
tection in a map may be employed to examine cluster selection (Bartlett 2000). 
It supposes that one can identify and then integrate over all the pixels covering a 
cluster (out to the virial radius) to access the total cluster signal Y for detection; 
it is therefore the best one can hope to do. 

Using this approach, the effects of resolution on Mjet are presented in Figure 
1, where the dashed lines (for the SZ curves) give M(jet(-z) once cluster extension 
has been taken into account. Clearly, clusters at low/intermediate redshifts are 
resolved out relative to the point -source calculation, because they are covered by 
several noisy pixels; this is most evident for the higher resolution interferometer- 
like curve. The difference for a low resolution survey, such as the Planck survey, 
is relatively minor, as can been seen. 

This point is important when it comes to predicting the number of objects 
to be expected. Clearly the predicted number will be lower when correctly 
accounting for resolution effects. Even more relevant to the question at hand 
(deep or wide), is the fact that, as seen in Figure 2 (and in Figure 4 below), 
the actual counts as a function of sensitivity are steeper than the point-source 



^For resolved sources, these counts are not necessarily the same as the counts as a function of 
integrated cluster signal Y . 
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Figure 2. Integrated SZ counts as a function of pixel noise expressed 
in terms of Y (arcmin^). Two cosmological models are used to illustrate 
the effects of source extension. The upper (red) curves are for an 
open model and the lower (blue) curves are for a critical model. The 
solid curves take into account finite source extension with a simple 
analytical model (Bartlett 2000), while the dashed lines give the results 
of the point-source approximation. In all cases, two different low-mass 

cutoffs - 10^3 ^ iq14 gQjg^^ 

masses - are applied in order to highlight 
the contribution of low-mass objects. The effect of the latter cutoff is 
clearly seen as the point where the counts flatten off at low Y. 
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calculation would suggest. In fact, at the bright end they are steeper than i^u~igg, 
the critical slope for answering our question. Counts steeper than —2 imply that 
more objects will be found by integrating deep on a small patch of sky (as long 
as l^oise ~ t'^^"^, which usually applies). At the faint end, the counts flatten 
out and rejoin the point -source calculation, which is always flatter than Y^f^^. 
Two important conclusions follow: 

1 . the point -source calculation gives the wrong answer in suggesting that one 
should go wide; 

2. the correct answer is found by accounting for resolution effects and indi- 
cates that, for each experimental set-up, there is an optimal sensitivity 
corresponding to the point where the counts reach a slope of —2. 



4. Expectations 

A number of dedicated instruments optimized for SZ surveying have recently 
been proposed. At the time of writing, many of these projects are funded and 
under construction, including: 

• Interferometer arrays 

— Arcminute MicroKelvin Imager (AMI) 

(tittp : //www.mrao . cam. ac.uk/telescopes/ami/ index .html) 

— Array for Microwave Background Anisotropy (AMiBA) 
( tittp : //www, asiaa. sinica. edu. tw/amiba/| ) 

— SZ Effect Imaging Array 

( |http : //www, astro .uiuc . edu'^jmohr/SZE/ index.html ) 

• Bolometer cameras 

— ACBAR 

( |http : / / cfpa.berkeley . edu'-^swlh/research/acbar .html ) 

— BOLOCAM 

( |http : //www, astro . caltech. edu^lgg/bolocam/bolocam.html ) 

• Planck satellite 

( |http : / / astro . estec. esa . nl/Planck/ ) 

Several authors have recently made predictions for these various instruments 
(Aghanim et al. 1997; Holder et al. 2000; Lo et al. 2000; Benson et al. 2001; 
Kay et al. 2001; Kneissl et al. 2001). In Figures 3 & 4 I show some expectations 
for the Planck Survey and for the AMI interferometer, both for two extreme 
cosmologies: a critical model and a purely open model with J7m = 0.3. A flat 
low-density model with non-zero cosmological constant will fall between these 
two cases, although somewhat closer to the purely open model. 

In the following, it should be kept in mind that the predicted total number 
of detected clusters is subject to large modeling uncertainties and can vary by as 
much as a factor of ~ 5. A particularly important source of uncertainty comes 
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from the uncertainty in the amplitude of the density perturbations erg. To be 
conservative, the predictions are given for 'low cs': {(Jg = 0.8, J^m = 0.3) and 
((Tg = 0.5, Om = 1) (e.g., Eke et al. 1998). Use of higher values quoted in the 
literature (e.g., Blanchard et al. 2000; PierpaoU, Scott, & White 2001) would 
substantially increase the predicted total number of clusters. Another source 
of uncertainty lies in the potential evolution of the cluster gas mass fraction 
/gas- The lack of evidence for strong evolution of /gas with redshift or for an 
important dependence on cluster mass may simply reflect a lack of sufficient 
data. In this light, it is important to note that the SZ signal for unresolved 
clusters (e.g., Planck detections) depends only on the global quantity /gas, and 
not on the details of the gas distribution; X-ray observations, as a contrast, are 
heavily influenced by the physics of the core, which is very poorly understood 
and difficult to model. 

The Planck predictions are given in terms of redshift distributions: shown as 
a function of redshift z are the total number of clusters with Y > lO"'' armin^ at 
redshifts larger than 2]^. The numbers at z = therefore give the total number of 
detectable clusters, and we see that on the order of 1 cluster /sq.deg. is expected 
for the survey, depending on the cosmology. At z ~ 1 the two models differ 
wildly. There are essentially no clusters expected in the critical model. On the 
other hand, about 10% of clusters are expected to lie at 2: > 1 in the open 
scenario. These conclusions assume that /gas remains relatively constant with 
redshift and mass. The presence of an important population of clusters at high 
redshift would be a strong indication of a low-density universe; but the lack 
of such a population could be accommodated in the open model by decreasing 
/gas with z. The result in any case would be an important piece of cosmological 
information. 

The predictions for an AMI-like instrument (resolution of a few arcminutes 
or better) are given in terms of integrated source counts as a function beam 
sensitivity. Both cosmological models are once again shown, each one for two 
different low-mass cutoffs (10^^ and 10^^ solar masses) in order to separate 
the contribution of low and high mass objects. The dashed curves show the 
equivalent counts in a point-source approximation, and we observe that the 
latter fails rather badly because most clusters are resolved at this resolution. 
The dotted lines in the background follow the critical slope of Y~'^. This kind 
of instrument reaches its optimal sensitivity at ~ few x 10~^ arcmin^/beam, 
where the predicted counts attain a slope of —2. At this sensitivity, one may 
expect, as a conservative round number and within the modeling uncertainties 
mentioned, on the order of several clusters/sq.deg. 

5. Conclusion 

It will shortly be possible to perform purely SZ-based surveys, opening up what 
may be refered to as the 'SZ band' for cluster studies. A number of ground- 
based instruments will realize such surveys within a few years over patches of 



■^i.e., the integrated counts from z — 00 down to z, as a function of 2 and at the fixed sensitivity 
of 10"'* arcmin^. 
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Figure 3. Planck-like cumulative redshift distributions for two ex- 
treme cosmologies. The y-axis gives the counts at y > 10"'^ arcmin^ 
for objects at redshifts > z, plotted as a function of z (x-axis). The up- 
per (red) curve corresponds to the open model, while the lower (blue) 
curve is for the critical model. In addition to the various modeling 
uncertainties mentioned in the text, the predicted number of clusters 
depends on the final sensitivity of the survey, which for the fixed ob- 
serving time of the mission will depend on various factors, including the 
weighting given to each of the bands, the effects of foreground removal, 
etc... This is an area of current study. 
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SZ Source Counts 




Figure 4. Integrated counts for an AMI-like instrument as a function 
of survey sensitivity in terms of noise rms/beam. The color code is the 
same as above. The thick sohd hues show the total counts, while the 
thin solid lines show only unresolved sources. Each model is plotted 
for two different low-mass cutoffs (10^^ and 10^^ solar masses); the 
higher mass cutoff in each case can be seen to flatten off at the low- 
noise end. The dashed lines show the results of the (invalid) point- 
source approximation, for both mass cutoffs and for the open model 
(above) and the critical model (below). Finally, the dotted lines in the 
background indicate the critical slope of —2, steeper than which deep 
surveys yield more objects than wide, shallow ones. 
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sky covering several square degrees; and by the end of the decade, the Planck 
mission will produce a catalog of several 10^ SZ clusters over the majority of 
the sky. Although predictions for the number of clusters expected for these 
surveys suffer from modeling uncertainties, various studies indicate roughly on 
the order of ~ 1 cluster /sq.deg. for Planck and ~ few clusters/sq.deg. for the 
deeper ground-based efforts under way. 

These surveys represent a new method for studying clusters that will find 
its place alongside X-ray and optical methods. As a population, clusters provide 
a unique statistical sample for probing cosmology and the large-scale structure 
formation process. Each observational approach to surveying for clusters has its 
proper advantages and disadvantages (collectively, its biases). In light of cosmo- 
logical studies, it is most instructive to quantify observational biases in terms of 
cluster mass and redshift, the two primary theoretical cluster descriptors. The 
specific advantages inherent to SZ surveying are listed above under Section 2 
and include the fact that the integrated SZ signal is proportional to a physical 
quantity - the ICM thermal energy. 

For a more complete understanding of the cluster population in view of its 
application as a cosmological probe it is essential to combine a variety of survey- 
ing (catalog selection) methods with different biases. SZ surveying will soon add 
an important new, independent and complementary approach in this endeavor. 

I would like to thank the organizers for their invitation to a most interesting 
and enjoyable meeting. 
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